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ABSTRACT 



Evidence is presented that indicates that Monterey Sub- 
marine Canyon was once the terminus of a major land drainage 
system. This pre-existing drainage system is not in evidence 
today because it has been altered by displacement along the 
San Andreas Fault. A mamerical model based on conservation 
of mass and plate tectonic reconstructions is utilized to 
reconstruct the topography of the region as it appeared prior 
to onset of motion along the San Andreas Fault. Model results 
indicate that the Colorado River may have drained into 
Monterey Bay during early Miocene time. 
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INTRODUCTION 



I . 

Monterey Submarine Canyon is the largest of the California 

submarine canyons. Traced seaward from its head near Elkhorn 

Slough in Monterey Bay, the canyon axis follows a meandering 

path for approximately 100 km before emerging as a fan valley 

at a depth of about 3 km (Figure 1) . From this point, the 

fan valley can be traced for another 400 km across the Monterey 

deep-sea fan. The volumes of the canyon and the fan have been 

3 3 

estimated to be 450 km and 30,000 km respectively (Menard, 
1960) . The walls of the canyon are approximately 1,500 m high 
at one point, producing relief comparable to that of the Grand 
Canyon of the Colorado River (Figure 2) . 

Such is the magnitude of this feature that two of its 
tributaries have been named. Carmel Canyon branches south 
and heads very close to shore off Monastery Beach in Carmel 
Bay. Sequel Canyon branches northward near the canyon head 
and extends toward the town of Sequel, east of Santa Cruz. 

The major canyon axis follows a fault contact near the 
head and is flanked by Tertiary sediments on both sides. At 
its outer limits; beyond about 20 km from its head at Elkhorn 
Slough, there is no indication of faulting along the canyon 
axis and the canyon course is purely erosional (Greene, 1977). 
Sequel Canyon is cut entirely in poorly indurated seimentary 
strata, while Carmel Canyon is eroded in granodiorite on the 
east and highly indurated sedimentary rocks on the west. 
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Figure 1. The Monterey and Carmel Submarine Canyons 

off the central California Coast. (Diagram 
by Tau Rho Alpha, USGS) 
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Figure 2. A comparison of the profiles of the Grand 

Canyon of the Colorado River and the Monterey 
Submarine Canyon, showing them to have similar 
relief. Elevations relative to sea level are 
given in feet. Vertical exaggeration is 5X. 
(After F. P. Shepard, Si±)marine Geology , 3rd 
Ed., Harper and Row, 1973.) 
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Like many of the submarine canyons along the California 
coast, Monterey Canyon does not lie seaward of a large land 
drainage system. The Salinas River is the largest of 
several small rivers which drain into Monterey Bay, all of 
which appear diminutive in comparison with Monterey Canyon. 

An explanation of this enigma was proposed by Martin and Emery 
(1967) when they suggested that Monterey Canyon had received 
drainage from the Great Valley of California, through the San 
Francisco Bay region during late Miocene, Pliocene and early 
Pleistocene time. Although some questions were answered by 
this hypothesis, others remain unanswered. Martin and Emery 
noted that the amount of material in Monterey Fan could not 
have resulted from the Great Valley connection which they 
proposed. In addition, they noted the existence of buried 
erosional features which led them to conclude that Monterey 
Canyon is a re-excavation of a pre-existing submarine canyon. 

As Martin and Emery noted, these buried erosional features, 
Elkhorn Erosion Surface and Pajaro Gorge, pre-dated their 
Great Valley Connection. The age and size of this pre-existing 
canyon were established by Greene (1977) using seismic profiling 
techniques. Greene's profiles confirm the fact that this pre- 
existing canyon existed prior to early Miocene time and is 
still not entirely re-excavated in some areas (Figure 3) . 

Greene ' s evidence indicates that the lower reaches of Monterey 
Canyon have been displaced northward and are represented today 
by Pioneer and Ascension Canyons. According to Greene's theory. 
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Figure 3. 



Basement contour map Monterey Bay Region, 
California (from Greene, 1977) 
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displacements along the Palo Colorado-San Gregorio Fault and 
the Ascension Fault over the past 20 million years may have 
moved these canyons into their present positions. 

Although it is presently one of the largest submarine 
canyons in the world, there is evidence to indicate that 
Monterey Canyon is an incomplete re-excavation of a larger 
canyon. The mechanism involved in excavation of this earlier 
canyon and creation of Monterey Fan is unknown. 

One possible explanation is that Monterey Canyon was the 
terminus of a large land drainage system which existed in 
pre-ear ly Miocene time. Large scale deformation of topography 
since early Miocene has removed the canyon from its source 
and produced changes in the drainage patterns of the area 
(Clark and Rietman, 1973; Greene, 1977; Martin and Emery, 

1967; Starke and Howard, 1968) . 

The topographic changes which have occurred since early 
Miocene along the west coast of North America have been exten- 
sive. About 30 million years ago, dioring Oligocene time, the 
East Pacific Rise came in contact with the North American 
Plate. Relative motion between the Pacific and North Ameri- 
can plates came to be expressed along right lateral strike- 
slip faults on the continental margin. During early Miocene 
time, the zone of strike-slip faulting between the Pacific 
and North American Plates shifted inland to the San Andreas 
Fault. Since its inception, approximately 300 km of right 
slip has occurred along the San Andreas Fault (Crowell, 1962) . 
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Such large-scale motion would undoubtedly have had great im- 
pact on any pre-existing drainage system. Reconstruction of 
topography as it existed prior to the onset of motion on the 
San Andreas Fault is difficult. Motion along the fault has 
not been constant in either speed or azimuth over this inter- 
val. In addition, compressional forces along the fault have 
uplifted areas resulting in the creation of new topographic 
features (Greene, 1977) . 

Due to the complexity of the problem, a mxmerical model 
is utilized to reconstruct the topography of the area as it 
appeared in early Miocene time. The model simply reverses 
motion along the San Andreas Fault utilizing data from plate 
tectonic reconstructions and imposing mass conservation. It 
does not reverse the effects of erosion or deposition, motion 
along other faults, or any other forces. However, since motion 
along the San Andreas Fault has contributed greatly to the 
alteration of topography since the early Miocene, results from 
such a model may be expected to provide useful insights into 
early Miocene topography . 
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II. PROCEDURE 



The model was designed to simulate a reversal of the 
motion which has occurred along the San Andreas Fault over 
the past 21 million years. In simulating this reversal, 
mass must be conserved. The conservation of mass states 
that : 



din 3m, ,,, 

g:^ = ^+ v*Vm + mV*v = 0, (1) 

where 

m = mass, 
t = time , and 
V = velocity. 

In other words, the total change in mass is equal to the local 
time change plus mass flux divergence. Since the total change 
is zero, mass is conserved and the local change must balance 
transport ; 

= - V • Vm - mV • V . ( 2 ) 

O U 

Since total mass equals mass density (p) times total volume 
(V) , then: 

= - V • V(pV) - pVV -V (3) 
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If we assume that density is constant in space and time, 
it can be eliminated: 

||.= _v-VV-VV-v (4) 

Finally, when dealing with a unit area, volume is ex- 
pressed in height above/below sea level (h) (Equation (5) . 

1^ = - V • Vh - hV • V (5) 

From Equation (5) , with an initial height field and a 
velocity field, topographic changes in time may be predicted. 

The height field is a 70 x 64 array of topographic and 
bathymetric heights as illustrated in Figure 4 , The grid 
interval is 10 minutes of longitude in the x direction and 
10 minutes of latitude in the y direction (15.2 and 18.5 km, 
respectively) . This interval in the x direction corresponds 
to true grid spacing at 35 degrees north, and introduces a 
departure from a spherical earth of approximately 7% at the 
northern and southern extremes. 

While the height field simply provides initialization for 
the model, the velocity field defines the fault and drives 
the model. As illustrated in the initial velocity field 
(Figure 5) , the portion of the grid which lies to the east 
of the San Andreas Fault is held fixed while that portion to 
the west of the fault is moved southeastward. Along the fault 
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Figure 4. Grid features 
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Figure 5. Velocity field: 0-4.5 million years before 

present, vectors indicate direction and 
relative magnitude of velocity at each point, 
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Figure 6. Velocity field: 4.5-10 million years before 

present, vectors indicate direction and 
relative magnitude of velocity at each point. 
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Figure 7 



Velocity field: 10-21 million years before 

present ; vectors indicate direction and 
relative magnitude of velocity at each point 



19 



itself, there is a gradation of the velocity field as depicted 
in Figure 8. At the position of the fault, the magnitude of 
the velocity field is reduced by one-half. At the next grid 
point to the east, the velocity is one-quarter as large as 
it is in the moving block. At the first grid point to the 
west of the fault, velocity is three-quarters as large as 
the velocity of the moving block. This configuration results 
in a zone of velocity gradation across the fault and prevents 
the numerical instabilities which may arise from modelling 
the fault as a more severe velocity discontinuity. Previous 
attempts without a gradation of velocity across the fault 
resulted in unreasonable topography. The components of the 
velocity field, speed and direction, will be discussed 
separately . 

Directions were derived from an algorithm which generates 
smooth fields to match directions at three points where direc- 
tions are known from plate tectonic reconstructions (Silver 
et al . , unpub. data) . Directions of relative motion between 
the North American and Pacific Plates relative to three points 
- on the North American Plate during three time intervals are 
presented in Table I. Corresponding values from the model 
algorithm are displayed in parentheses for comparison. Since 
the model simulates a reversal of motion, all directions are 
reduced by 180° to produce the velocity fields depicted in 
Figures 5-7. These figures illustrate the velocity fields 
used at different times in the model. The vectors to the left 
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TABLE I 

(From Blake al . , 1978) 

Azimuth of Average Movement Direction of Pacific 
Plate Relative to Points on North American Plate 

Values from model algorithms in parenthesis 



Age (my) 


36°N, 121. 5°W 


33°N, 119°W 


26°N, 


112°W 


0 — 4.5 


321.2° (321.2°) 


318.6° (318.5°) 


311.7° 


(311.7°) 


4.5—10 


328.0° (328.0°) 


325.7° (325.6°) 


319.5° 


(319.5°) 


10 — 21.2 


339.0° (339.2°) 


335.0° (334.9°) 


323.9° 


(324.0°) 



side of the figures represent the moving block. It can be 
seen that there is a very slight change in the direction from 
one end of the field to the other in all three figures as 
would be expected from Table I. On the right side of each 
figure is an area with no vectors. This is the stationary 
portion of the grid. The shorter vectors at the intersection 
of the stationary block and the moving block represents the 
fault zone. As previously discussed, there is a gradient in 
the magnitudes of the velocities across the fault. This is 
reflected in the comparatively shorter vectors in the fault 
zone. The actual position of the San Andreas Fault is super- 
imposed on these figures for comparison. It can be seen from 
Figures 5-7 and from the values in Table I that the motion 
changes with time to a more westerly direction becoming less 
compressional and more nearly aligned with the fault axis. 

Speeds were derived from a combination of sources. Plate 
tectonic reconstruction (Atwater, 1973) produces an accelerating 
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‘ velocity field which would have resulted in approximately 
I 610 km of offset if all motion had been expressed as fault 
displacement. How much of this displacement is actually 
expressed in offsets along the San Andreas and other faults 
is unclear. Evidence for a 260 km offset along the San Andreas 
Fault (Crowell, 1962) establishes a conservative estimate 
and is the reference displacement utilized in this model. Al- 
though displacements from plate tectonic reconstruction do not 
match displacements derived from geologic evidence, the accelera- 
tion of motion between two plates, as derived from plate tec- 
tonic reconstruction, is not disputed by geologic evidence 
(Atwater, 1973). This acceleration of motion is incorporated 
into the model by reducing velocities from plate tectonic 
reconstruction such that the total resultant offset equals 260 
km. This is accomplished by multiplying the velocities by 
a scaling factor, s = 260 km/610 km = .426. Resultant scaled 
velocities as well as original velocities derived from plate 
tectonic reconstruction (Atwater, 1973) are listed in Table II. 



TABLE II 



Time Period 



Relative Motion 



Relative Motion 



(Million years 
before present) 



Pacific/N. Ameri- 
can Plates 



San Andreas 
Fault 



21—10 



1 . 3 cm/yr 



0 . 6 cm/yr 



10—4.5 



4 . 0 cm/yr 



1 . 7 cm/yr 



4.5— 



5 . 5 cm/yr 



2 . 3 cm/yr 
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The height and velocity fields are substituted into the 
Continuity Equation (Equation (5)) using finite differencing 



schemes. In discussing differencing schemes, the following 
coordinate system and symbology will be used: 



H, . . , . 

(i,D /t) 



u 



reference point, height at a location on 
a horizontal plane identified by coor- 
dinates i (east) and j (north) at time t 

velocity component in the i (east) 
direction 



V = velocity component in the j (north) 
direction 



At = interval of time (200,000 years in this 
model) 



Ax 



Ay 



interval of space on the i axis . 
example, the distance between H 



For 
(i, j /t) 



interval of space on the j axis. The 

distance between H,. . ,, and H,. 

(i/3,t) (i,]+l,t) 



The primary centered differencing scheme used in this 
model expresses Equation (5) in the following way: 



H(i, j,t+l) - H(i,i ,t-l) 
2At 



-u(i, j ,t) [H(i+1, j , t) -H (i-1, j , t) ] 

2Ax 

v(i,j,t) [H (i , j + 1, t) -H (i, j-1, t) ] 

2Ay 



H(i,j,t) [u(i+l, j ,t)-u(i-l, j ,t) ] 

2Ax 



H(i,j,t) [v(i, j+l,t) -v(i, j-1, t) ] 

2Ay 



( 6 ) 
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This scheme is centered in space and time. It can not be 
used initially since it requires height fields at two differ- 
ent times (t and t-1) in order to predict the height at t+1 . 
In addition, it can't be used on the grid boundaries since 
it would require the input of values beyond the range of the 
data field. However, it offers the advantage of centered 
differencing and favorable numerical stability. 

Forward differencing must be used to predict the height 
fields on the first iteration for sxibsequent input into the 
primary centered scheme. Using forward time and centered 
space differencing. 



H(i,i,t-H) -H(i, j ,t) 
At 



u(i,j,t) [H(i-H, j ,t) -H(i-1, j ,t) ] 
2Ax 



v(i,j,t) [H (i, j-t-1 , t) -H (i, j-1, t) ] 
2Ay 



H(i,j,t) [u(i+l, j , t) -u(i-l, j ,t) ] 
2Ax 



H(i,j,t) [v(i, j + 1 , t) -v(i, j-1, t) ] 
2Ay 



(7) 



As with all mamerical methods, numerical stability is an 
important consideration and limits the size of the time step 
(At) to be utilized. The Courant-Friedrichs-Levy (CFL) 
condition (Courant ^. , 1928) for computational stability 
(Equation (8)) applies for this model. 



V At 
Ax 



< 



1 



( 8 ) 
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It controls the size of the time step (At) and grid 
spacing (Ax, Ay) such that motion does not cover more than one 
grid space between computations. With this model, At must be 
less than or equal to 500,000 years. The time step used is 
200,000 years. 

Boundary conditions are illustrated in Figure 4. As pre- 
viously discussed, centered space differencing can not be 
used on the perimeters of the grid. Two different types of 
boundary conditions are used in this model. Boundaries in the 
direction of motion of the moving block (the eastern and southern 
boundaries) must be permeable. For this reason, the Upstream 
Differencing scheme (Haltiner and Williams ,_ 1980 , p. 130) is 
used on these boundaries. This scheme produces a permeable 
boundary where it intersects the moving block and establishes 
a rigid boundary condition for the stationary block: 



H (i, j , t+1) -H (i, j , t) 
At 



u (i-1, i , t) [H (i,i , t) -H (i-1 , i , t) ] 
AX 



H(i-l,j,t) [u(i, j ,t) -u(i-l, j ,t) ] 
Ax 



( 8 ) 



Fixed boundary conditions are imposed on the northern 
and western boundaries . Maintaining the points on these 
boundaries at their initial heights introduces errors where 
these points lie on the moving portion of the grid. As 
topography on these boiandaries moves toward the interior 
of the grid, a false field is created in their previous 
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positions. Resultant erroneous values are eliminated from 
output by moving the northern and western boundaries inward 
(for display only) such that the northwest corner of the grid 
maintains its position relative to the moving block. This 
is illustrated in Figure 9 . 

The model produces output at 3 million year intervals. 
Output is displayed in a three dimensional form as well as a 
contoured plan view at each interval (Figures 10-25) . The 
three dimensional computer graphics simulate viewing the ter- 
rain from a vantage point at high elevation over the Pacific 
Ocean looking northeast. Submarine bathymetry is included and 
all vertical heights are exaggerated by a factor of 15. 
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Figure 9. Movement of boundaries for display 
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Figure 10. Topography 21 million years before present 
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Figure 11. Twenty-one million years before present. On 
all contour plots; contour interval = 500 m, 
figures in hundreds of meters, bathymetry in 
dashed lines, inferred drainage is heavy 
dashed line, M, B, and L represent Monterey, 
Bakersfield and Los Angeles respectively 
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Figure 12. Topography 18 million years before present 
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Figure 13. Eighteen million years before present, 
contour interval = 500 m, figures in 
hundreds of meters 
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Figure 14. Topography 15 million years before present 
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Figure 15. Fifteen million years before present, 
contour interval = 500 m, figures in 
hundreds of meters 
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Figure 16 . Topography 12 million years before present 
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Figure 17. Twelve million years be;fore present, contour 
interval = 500 m, figures in hundreds of 
meters 
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Figure 18 . Topography 9 million years before present 



37 




Figure 19. 



Nine million years before present, contour 
interval = 500 m, figures in hundreds of 
meters 
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Figure 20. Topography 6 million years before present 
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Figure 21. Six million years before present, contour 
interval = 500 m, figures in hundreds 
of meters 
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Figure 22. Topography 3 million years before present 
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Figure 23. Three million years before present, contour 
interval = 500 m, figures in hundreds of 
meters 
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Figure 24. Present day topography 
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Figure 25. 



Present day, contour interval 
figures in hundreds of meters 



500 m, 
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III. RESULTS 



In viewing the figures depicting model reconstructions 
of paleotopography , the design of the model must be kept in 
mind. Coarse resolution, as well as the simplifying assump- 
tions of the model limit the scope of consideration to large 
scale effects. 

Figures 10-25 are model outputs generated at three 
million year intervals. At each interval there is a three- 

k 

dimensional topographic display and a contour map, both 
constructed from model results for that particular time. 

In all three-dimensional plots, the coastline (based on 
present sea level) is represented by a heavy solid line. 

The points corresponding to the cities of Monterey, Bakers- 
field, and Los Angeles have been annotated for reference. 

All contour maps are oriented with north at the top. 
Bathymetry appears as dashed contours while all contours ^0 
appear as solid lines. The San Andreas Fault is represented 
as a heavy solid line. Inferred drainage patterns are 
indicated by a heavy dashed line. The points corresponding 
to the cities of Monterey, Bakersfield, and Los Angeles are 
labelled M, B, and L respectively. 

Although the model runs backward in time, this discus- 
sion will start with model output at 21 million years before 
present and proceed foward in time for simplicity. As 
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output for each time interval is discussed, reference will 
be made to both figures (3-D and contour) representing 
output for that interval . 

As expected, topography at 21 million years before 
present (Figures 10 and 11) differs greatly from present- 
day topography. Monterey Bay lies far to the south-east 
of its present position, approximately 100 km west of 
Bakersfield. Bakersfield, and points to the east of the 
San Andreas Fault, appear much as they are today. The Great 
Valley of California extends northwest of Bakersfield and is 
bounded on the east by the Sierra Nevada Mountains. The 
Mojave Desert lies to the south of the Sierra Nevadas . 

Along the fault itself, topography is depressed, particularly 
in the area to the southwest of the Mojave Desert. Areas to 
the west of the fault appear radically different than they 
do today. All topographic features which correlate to today's 
topography appear far to the southeast of their present posi- 
tions. Although the general shape of the coastline and the 
mountain ranges are recognizable, they appear different than 
they do today. Monterey and San Francisco Bays are different 
in shape and location but are easily recognizable. The Coast 
Ranges to the north appear at lower elevations than they do 
today but are not greatly different in the area to the south 
of Monterey. The Transverse Ranges are somewhat lower and 
they are located to the southeast of their present positions 
as are the Peninsular Ranges. The position of the Peninsular 
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Ranges correlates geographically with the present-day 
position of the Gulf of California. 

Despite the fact that this reconstruction is the product 
of a relatively simple geological model which does not take 
into account the effects of erosion, deposition, or displace- 
ments along other faults, the large scale features which 
appear in these figures are supported by geologic evidence. 
The Sierra Nevada Mountains (Curtis et a^. , 1958) and the 
Colorado Plateau (Eardley, 1962) are features which came 
into existence long before Miocene time. Also, the Gulf of 
California did not exist at this time (Larson ^ , 1968; 

Larson, 1971; Van Andel and Shor, 1964). The depressed 
areas along the fault southwest of the Mojave Desert corre- 
late geographically with the Salinas, Caliente, San Joacquin, 
Ridge, and Soledad basins. The geologic history of these 
basins (Norris and V7ebb, 1976; Blake et ad., 1978) also 
shows good chronological correlation with the model 
reconstruction . 

It is apparent from Figures 10 and 11 that drainage 
patterns were much different in early Miocene time. Drainage 
from the Colorado Plateau (Colorado River) could not flow 
southward into the Gulf of California as it does today. 

Rather this drainage flowed westward and entered the ocean 
somewhere to the north of the Transverse Ranges. Exactly 
where the terminus of this drainage system was located can 
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not be determined from model results alone. However, evidence 
previously presented indicates that Monterey Bay was the 
terminus of a major land drainage system at this time (Greene, 
1977). The inferred drainage is illustrated in Figure 11. 

How long this drainage pattern had been in existence prior 
to early Miocene time is a question beyond the scope of this 
study. However, insights into changes which have occurred 
since early Miocene may be obtained by examining model output 
at intervals over the last 21 million years. 

As expected from the relatively low velocity of motion 
along the fault from 21 million years before present to 10 
million years before present (Table II) , there is little 
change in topography over this period (Figures 10-17) . It 
appears that the drainage pattern previously established 
would have been preserved over this interval. Although there 
is evidence to indicate that Monterey Submarine Canyon was 
filled and re-excavated during this time (Greene, 1977) , this 
could have been the result of sea level fluctuations and does 
not necessarily imply large scale changes in land drainage 
patterns . 

From 10-4.5 million years before present, motion along 
the fault increased in speed and became slightly less com- 
pressional in azimuth (Tables I and II, Figure 6). Topography 
at 9 million years before present (Figures 18 and 19) shows 
that compression along the fault in the area southwest of the 
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Mojave Desert is closing the depressions in that area. At 
the same time, the Peninsular Ranges have been moving north- 
westward, leaving a depressed area in the southeastern corner 
of the grid. At this point it is not possible to determine 
drainage patterns, as illustrated in Figure 19. At 6 million 
years before present (Figures 20 and 21) , it is clear that 
compressional motion along the fault has uplifted the area 
to the southwest of the Mojave Desert to such an extent that 
this area no longer serves as a conduit for drainage. Simul- 
taneously, the depressed area in the southeast corner has 
expanded due to continued movement of the Peninsular Ranges 
out of this area. It is not possible to determine that 
drainage actually flowed southward at this point since 
topography beyond the southeastern corner of the grid is not 
presented. However, geologic evidence (Larson et ^. , 1968; 
Larson, 1971) indicates that the Gulf of California developed 
at about this time. For this reason, it is assumed that 
drainage from the Colorado Plateau flowed southward much as 
it does today. 

Velocity of relative motion is slightly higher over the 
last 4.5 million years to the present (Table II) and there 
is a slight change in the azimuth of direction (Table I and 
Figure 5). However, the trends concerning drainage continue 
during this interval as depicted in Figures 22-25. North- 
westward displacement of points to the west of the San 
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Andreas Fault, compression in the area of the Transverse 
Ranges, and expansion of the Gulf of California continue to 
the present and are continuing today. 
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IV. CONCLUSIONS 



In arriving at conclusions based on the results of a 
numerical model, care must be taken to consider the inherent 
weaknesses of the model and to examine available evidence 
from other sources for verification. In presenting model 
results, only large scale features and general trends were 
discussed. It is clear, from the description of the model, 
that these are the limits of its credibility. However, 
within these limits, model results correlate well with 
available evidence and indicate that the Colorado River was 
the erosional force involved in the excavation of Monterey 
Canyon . 
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